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Abstract The reorientation of the magnetization of a ferromagnetic monolayer is calculated with the help 
of many-body Green's function theory. This allows, in contrast to other spin wave theories, a satisfactory 
calculation of magnetic properties over the entire temperature range of interest since interactions between 
spin waves are taken into account. A Heisenberg Hamiltonian plus a second-order uniaxial single-ion 
anisotropy and an external magnetic field is treated by the Tyablikov (Random Phase Approximation: 
RPA) decoupling of the exchange interaction term and the Anderson- Callen decoupling of the anisotropy 
term. The orientation of the magnetization is determined by the spin components {S°') {a = x, y, z), which 
are calculated with the help of the spectral theorem. The knowledge of the orientation angle &o allows 
a non-perturbative determination of the temperature dependence of the effective second-order anisotropy 
coefficient. Results for the Green's function theory are compared with those obtained with mean-field 
theory (MFT). We find significant differences between these approaches. 
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1 Introduction 



, Experimental and theoretical investigations of the mag- 
' netic properties of ultrathin ferromagnetic films are a topic 
I of intense current interest pi. In particular, the temper- 
■ ature dependent (effective) magnetic anisotropy and the 
resulting direction of the magnetization have been deter- 
, mined for a variety of thin film systems |3]. The differ- 
' ent temperature and thickness dependence of the various 
anisotropy contributions may result in a reorientation of 
the magnetization as function of temperature and/or film 
' thickness. In addition, a magnetic reorientation can be 
, induced by applying a magnetic field @,Ql, which is one 
of the experimental methods for measuring the magnetic 
anisotropics. In order to determine the different contribu- 
tions to the magnetic anisotropy separately, their depen- 
dence on the temperature, on the film thickness, and on 
the magnetic field has to be known. 

By use of ab-initio methods, the magnetic anisotropy 
of thin films has been calculated as the total energy differ- 
ence between different directions of the magnetization g . 
Until now this has been done only for T = 0, whereas mea- 
surements are always perfornred at finite temperatures. 
Thus, the knowledge of the temperature dependence of 
the anisotropies is also required in order to allow a com- 
parison of experimental and theoretical results. 



This knowledge is gained mainly within the framework 
of a Heisenberg model for localized spins, which, how- 
ever, seems to yield satisfactory results also for the im- 
portant case of itinerant magnets such as Fe, Co, and Ni. 
Usually mean-field theory (MFT) is applied to calculate 
the magnetization, accompanied either by diagonalization 
of a single-particle Hamiltonian Iq], or by a thcrnrody- 
namic perturbation theory for the anisotropy terms [Q. 
In principle, MFT is not applicable to two-dimensional 
(2D) magnetic systems. In such systems, thermodynamic 
correlations, which are neglected by MFT, have a deci- 
sive influence on the magnetic properties. In particular, 
the long-range magnetic fluctuations destroy the rema- 
nent magnetization of an isotropic 2D Heisenbergmodel at 
finite temperatures (Mermin- Wagner theorem ||8|). How- 
ever, even small anisotropic contributions, which are al- 
ways present in real magnetic systems, induce a magneti- 
cally ordered state in thin films with a critical temperature 
of the order of the exchange coupling [|| . This is the rea- 
son why results obtained from MFT are expected to be 
qualitatively correct also for 2D Heisenberg magnets. Fur- 
thermore, it is known that applied magnetic fields have a 
much larger impact on the magnetic properties of 2D than 
of three-dinrensional (3D) systems |10|| . Therefore, mea- 
sured magnetic quantities such as the effective anisotropy 
may in principle depend on the experimental situation. 



p. Frobrich ei al: Field-induced magnetic reorientation and effective anisotropy 



e.g. whether a magnetic field is present or not. The MFT 
method poorly reproduces this sensitive field dependence 
of 2D systems, thus necessitating an improved theoretical 
description. 

Long-range magnetic fluctuations can be treated by 
spin wave theories. To our knowledge, the calculation of 
a field-induced magnetic reorientation with such an ap- 

§ach has been performed only by Erickson and Mills 
, who consider an exchange coupling, a uniaxial lat- 
tice anisotropy, and the magnetic dipole coupling. They 
have transformed the spin operators into Bose operators 
(Holstein- Primakoff transformation |Q), which are treat- 
ed only in lowest order. Thus the validity of this linearized 
spin wave theory is limited to low temperatures. The au- 
thors obtained a strong increase of the transverse fluctua- 
tions near the reorientation transition, where the direction 
of magnetization turns into the field direction. However, as 
remarked by the authors, the method should break down 
in this region. 

In a previous paper |l3[ , we demonstrated with the 
help of a many-body Green's function theory that the 
Tyablikov (or RPA) decoupling |lj] of the higher-order 
Green's functions provides a significantly improved de- 
scription of the magnetization over MFT. We showed this 
by comparing RPA and MFT results with the exact solu- 
tion of a Heisenberg spin pair, and with the 'exact' Quan- 
tum Monte Carlo result flq] of a Heisenberg monolayer 
with spin S = 1/2. Also the temperature dependence of 
the second- and fourth-order effective anisotropy coeffi- 
cients was calculated by a thermodynamic perturbation 
theory, expecting that RPA gives an improved descrip- 
tion also for the effective anisotropy coefficients. As re- 
sults we found that their temperature dependence looks 
different, particularly at low temperatures, and that their 
dependence on the magnitude of the spin is much weaker 
in RPA than in MFT [^. Such a perturbative approach 
makes sense only if the magnetic field is stronger than the 
anisotropy. 

In the present work, we investigate also in the frame- 
work of many-body Green's function theory the orienta- 
tion of a ferromagnetic monolayer at finite temperatures, 
the field-induced magnetic reorientation, and the effective 
(temperature-dependent) anisotropy. Encouraged by the 
fact that RPA yields a good approximation to the mag- 
netization, we expect that the RPA results of the present 
paper also represent more satisfactory estimates for the 
above mentioned magnetic properties than the results of 
MFT. Added to the Heisenberg Hamiltonian is the second- 
order single-ion anisotropy favoring a perpendicular mag- 
netization, and an external magnetic field perpendicular 
to this uniaxial lattice anisotropy, causing a magnetic re- 
orientation with an increasing field strength. Since we are 
mainly interested in the action of the second order single- 
ion anisotropy and its temperature dependence, we omit 
here the magnetic dipole coupling. For a ferromagnetic 
monolayer the dipole coupling competing with the out- 
of-plane second order single-ion anisotropy will induce a 
reorientation as a function of the temperature only in a 
narrow parameter range (for the strengths of the dipole 



coupling and the anisotropy) and therefore a reorientation 
is quite improbable in general, see e.g. |lq| . The magne- 
tization of a monolayer will in most cases stay either in- 
plane or out-of plane. This is the reason why we study the 
reorientation induced by an external magnetic field. For 
several layers, on the other hand, the dipole coupling will 
play a more important role. A magnetic reorientation as 
function of film thickness and/or temperature can be ex- 
pected, when the film surface and film interior anisotropic 
contributions compete with each other and exhibit differ- 
ent temperature dependences. 

We apply the RPA method for decoupling the Green's 
functions coming from the exchange coupling terms thus 
approximately taking into account interactions between 
magnons, whereas we make use of the Anderson-Callen 
decoupling p7| for the corresponding anisotropy terms. 
The magnetization axis will be tilted with respect to the 
easy axis by the magnetic field; this procedure resembles 
the experimental situation [||g] . We do not rotate the lo- 
cal spin quantization axis but calculate the magnetiza- 
tion from the expectation values of the spin components 
(5'"). This gives the equilibrium orientation angle directly. 
Knowledge of the orientation angle allows a calculation of 
the temperature dependent (effective) anisotropy coeffi- 
cient from the condition that the free energy be a min- 
imum. This is a non-perturbative approach because the 
quantities entering in the final expression are calculated 
from the full Hamiltonian. This is an improvement over 
the previously ||I^ used thermodynamic perturbation the- 
ory, where the anisotropy term is treated as a small per- 
turbation. We compare the effective anisotropy obtained 
from the non-perturbative approach with that from the 
thermodynamic perturbation theory and, in many places 
throughout the paper, we compare the results obtained 
from the Green's function theory with those from MFT. 



2 The Green's function formalism 

Our aim is to determine the orientation angle 6*0(2^) of 
the magnetization of a single (001)- layer as a function 
of the temperature T from the expectation values (S"") 
(a = X, y, z) of the components of the magnetization. For 
this purpose a Heisenberg Hamiltonian is used [^|7| con- 
sisting of the isotropic exchange interaction J between 
nearest neighbour lattice sites, a second-order single-ion 
anisotropy parameter K2 = K2(T — 0) at zero tempera- 
ture, and an external magnetic field, B = {B^ , B^ , B^) 



n 



J 



^j2^,s,-K^J2iSkr-E^s, 



kl 



J 



Y^is^s^ + s^sn^K^Y^isif 



kl 



J2(Ib-s+ + Ib+s^ + b^s{), 



(1) 



where the notation Sf ^ Sf ± iS^ and B^ = B"" ± iBV 
is introduced. 
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In order to treat the problem, one needs the foUowing 
Green's functions: 



terms do not vanish. Instead, an ansatz of the foUowing 
form is introduced |E3l 



G'r/(':r = ((^r;(^ir(^7n); 



(2) 



f: ±.?nn/-» lb, mn 



where ry = ±1 refer to the commutator {rj — —1) or anti- 
commutator [rj = 1) Green's functions, respectively, n > 1 
and m > Q are positive integers, i and j denote lattice 
sites. 

The G'^-T^T are determined from the equations of mo- 
tion in the spectral representation 

^ Gf;™"(c.) = A'^:^^-+{{[sf,nU; is^riSiT)). , (3) 

with the inhomogeneities 

= {sns;r{s-r + v{s]ris-rs^) , (4) 

where (...) = Ti{...e-^'^) with (3 = 1/fcsT and fcs Bohz- 
mann's constant. 

Knowledge of the Green's functions allows the determi- 
nation of the resp ective correlation functions by the spec- 
tral theorem |14| 






{{SrS! + 5f 5f ; iS^y-iSrr)) c <?f — G^'™", (8) 



where the functions <?j ' have to be determined. This 
type of decoupling is valid for anisotropics which arc small 
compared to the exchange interaction. In Appendix A, 
we have investigated various decoupling schemes for the 
single-ion anisotropy proposed in the literature. In the ap- 
pendix, we give arguments for our preferring the Anderson- 
Callen decoupling for the present calculations, yielding the 
decoupling function for n = 1 



<l> = (p 



±,?nl 



^{S!)[i^^[sis + i)-{s:sf)] 



Ctr = {{S'-riSjTST) = - Ihn 



dij 



cflul 



(9) 

We now apply the Tyablikov decoupling to the ex- 
change interaction term and the Anderson-Callen decou- 
pling to the single-ion anisotropy term in equations (ph. 
Performing in addition a 2D- Fourier transformation to 
momentum space, with k = (kx,ky,0) being the in-plane 
wave vector, yields the following set of equations of motion 

(uj-H^ i7+ \ /G+^""(k,w)~ 
u + H'- -H- G-^""(k,w) 
iff+ uj y yG^''""(k,u;) 



ii/- 



A -\- ,'17171 

A — .run 
Az.mn 
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ijiv) 






(10) 



(5) 



Calculation of the commutators [Sf ,-ff]_i yields the fol- 
lowing set of equations of motion for the Green's functions 



wG 



ij{r}) 
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{mst 



^k ^i 



1 



1 



{ QZ\m( Q—\n\\ ^ j:)— y-~i-\-,nui 



-b+g: 



(6) 



with the abbreviations 

ff"=B" + (5")J(q-7k), a = +,-,z 

H'- =H''+K2$ = Z + {S')Jiq - 7k) 
Z = B^ +K2^. (11) 

For a square lattice, one obtains 7k — 2(cosfc:r + cosky), 
and 5 = 4 is the number of nearest neighbours. Note that 
j^a.mn _ j^a,mnf]^-^ dcpcuds On the wavc vcctor k for 77 = 1 
but not for 77 = —1. 

The determinant of the matrix in equation (nfl) is given 

by 

A{uj,k)^iu{uj-Ei,){iu + Ei,), (12) 

which is the magnon dispersion relation 

Ek = Vh+H- + H^H^ . 



The higher-order Green's functions occuring on the right- 
hand sides have to be decoupled in order to obtain a closed 
set of equations. For the exchange coupling terms we ap- 
ply a generalized Tyablikov- (or RPA-) H,Q decoupling, 
allowing also for a finite value of the x, y- (or ±-) compo- 
nents of the magnetization (a, /? = -t-, — , z; i 7^ fc) 



Hence, the eigenvalues are 

ll>i =0 , a;2,3 = ±Ek . 
The Green's functions are given by 

'(u;,k) 



G° 



\i,,k) 



Aa^mnf 



A{LU,k) 



a = +,-,z, 



(13) 

(14) 

(15) 



{{SrS^,; (S^riSrr)) c. (5r)Gf:"™ + (5f )G 



~ia,mn 



(7) 



The terms resulting from the single-ion anisotropy (i — k) 
have to be decoupled differently. A RPA-decoupling, as 
was proposed by Narath |2^ , is reasonable for an exchange 
anisotropy (ex SfS^). It yields for the single- ion anisotropy 
unphysical results, for instance for S' = 1/2 the respective 



where Z\"'™" is the determinant where column a of Z\ 
is replaced by the right-hand side of equation (Kw. One 
obtains, for example 



i^'""(w,k) = A^^""^^ _ (Hz)2\ 



\ A+ '"'''' H-{uj + H-) ~ ^A--""'H+{iu - H-) . (16) 
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Now we use the fact that the commutator Green's func- 
tions must be regular for cj — > 0, e.g. ||l4|, 



Hm ujGt'^"'{uj, k) = . 



(17) 



Thus Z\"'™"(0,k) ~ 0. Since one of the eigenvalues van- 



ishes, see equation (|14|), we obtain from equation (16) for 



^-^+,m„ ^ H+AZf" + 2 iJM^f " = . 



(18) 



Evaluating this expression for m = and n — 1 we 
find together with the definitions in equation ([L1|) 






B^ 



Z 



B^ +K2<P 
Putting this into equation (p^ ) we have 



-2ZA' 



= A^f "B 



A" 



'B+ 



(19) 



(20) 



Equations (|lj) and (|20| ) are important relations between 
correlation functions which we call the regularity condi- 
tions. We note that the same relations are obtained for 
all three a = +,—, z. With m — and n = 1 we obtain a 
relation between (5'''} and (5^) 



(5±) 



B^ 



■is^) 



(21) 



B± {S') 



B' + 27^2(5^) (1 - [5(5 + 1) - (5^5^)]/252) 



This means that the knowledge of (5^) and {S^S^) also 
determines the expectation value of the spin components 
(5^) and (5^). Other useful relations obtained from equa- 
tion (20) are given in Appendix B. 

It remains to establish equations which determine the 
moments ((5^)™). For this purpose we consider for exam- 
ple the following commutator Green's function 



-1 i^, k) 



Z\l-r(c^,k) 
^(^,k) 



[uj - Ek){Lo + E^) 



LO 



A-^r"-H-uj + A-\""'-H+Lo 
2 ^1 2 



-1 



A 



z,7nn 2 
UJ 



(22) 



where the regularity condition has been taken into ac- 
count. 

In order to calculate correlation functions in the case of 
a vanishing eigenvalue one also needs the respective anti- 
commutator Green's function G^'™" , since in this case the 
correct form of the spectral theorem reads, e.g. plj 



— lim 

2tt s^a 



duj 



CO 

z,mn 



;/3^_l 



G'^T (^ +iS, k) - G!_f"(cc; - iS, k) 



Dt 



where 



Dt 



1™ t: G+'""(w,k). 



(23) 
(24) 



Using 



+1 (w,k) 



zA;f"(co,k) _ 1/2 

Aiuj, k) ^ uj{lu - Ei,){uj + Sk) 



^+,rnn H' {u + H^) + ^1;^ " H+ {uj - H^) 



2A;f" cj^-(iJ- 



(25) 



and the relation between anti-commutator and commuta- 
tor correlation functions 



A+i (k) = ^ Ji 



o /^amn 



(26) 



we find together with the regularity condition 

^k"" = ^(c+""'H-H^+C^""'H+H-+2C^"'"{H^')^ 
2-E^k V 

Finally, we obtain from equation ( P3D together with equa- 
tions (|2|) and {^ 



-A+'r'^EkH- 
2-1 

2-1 



-^ - coth(f3Eu/2) 
-^ + coth{f3Eu/2) 



(28) 



From these relations and the regularity conditions, equa- 
tion (|20|), one obtains all necessary expectation values. 

In the following we restrict ourselves to an external 
magnetic field B confined to the a;z-plane. Because of the 
azimuthal symmetry in the case of an uniaxial anisotropy 
it is sufficient to deal with the z- and x- components of 
the magnetization ((5^) = for B^ = 0). From these 
values the magnitude and the equilibrium polar angle of 
the magnetization are determined: 



M^T) = (5^)2 + (5^)2, 
Oo{T) = arctan = arctan 



gx 



(5^) 



B'- + K2<P 



(29) 



The knowledge of M{T) and 6*0 enables a non-pertur- 
bative determination of the temperature dependence of 
the anisotropy coefficient. The anisotropic part of the free 
energy is usually written as a power series of the direction 
cosines of the magnetization. Its precise form is written 
in accordance with the symmetry of the system and usu- 
ally converges. This series need not only be valid for small 
fields and anisotropics as compared to the exchange cou- 
pling. It can even be considered as a definition of the effec- 
tive (temperature-dependent) anisotropy coefficients. It is 
such an expression from which experimentalists determine 
the effective anisotropics. The corresponding part of the 
free energy in the presence of the lowest order term, the 
second-order uniaxial anisotropy K2{T), and the Zeeman 
term has the form 

F{T, O) = Fq{T) - K2{T) cos^ 6> - B ■ M(T) . (30) 
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This corresponds to the physical situation where the higher 
order anisotropics can be neglected. We only treat this 
case in the present paper. K2(T) is then determined from 
the condition that the free energy has a minimum at the 
equilibrium angle Oq: 

^^^ {0o)=O = 2K2 (T) cos 00 sm 0o 

- B^'MiT) cos 6>o + B^ Af (T) sin Oq . (31) 
Solving with respect to K2{T) we find 

""'^^^ = 2 cos 00 sin 00 ■ ^^'^ 

Note that the effective anisotropy can be considered as 
an intrinsic property of the layer itself only if K2{T) as 
obtained by equation (32) is practically independent of 
the external magnetic field. We mention that the effective 
anisotropy at zero temperature K2(T = 0) differs from the 
K2 in the Hamiltonian by a spin-dependent normalization 
factor §,|l|. By expanding equations (|), (|9|) and (H) 
for small T and 0o one finds K2iT = 0) = K2 S{S - 1/2). 
The results shown in figures 3, 4, 6 are normalized with 
respect to this factor for S = I. 

The procedure described above is non-perturbative in 
the sense that the quantities (the magnetization and the 
orientation angle) entering in equation ( p^ are calculated 
from the full Hamiltonian. This is in contrast to a thermo- 
dynamic perturbation theory in which the lattice anisotro- 
py is used as a small perturbation. The Hamiltonian is 
split into two terms Ti. = Ho — K2 J2i{^i)'^j ^^^ the effec- 
tive anisotropy is calculated from moments of the unper- 
turbed Hamiltonian Ho only, see e.g. Ref. (p3[). 

In the present thermodynamic approach K2{T) de- 
pends on the temperature mainly through the magnetiza- 
tion, which itself is a function of the magnetic field. Thus 
at least a weak field dependence of the effective anisotropy 
is expected. There are other sources which might induce 
a temperature dependence. In Ref. |T9] it has been shown 
within a tight-binding scheme that at a constant magneti- 
zation, the electronic entropy, the thermal lattice expan- 
sion, the population change of spin-orbit-splitted energy 
levels near the Fermi energy cause a decrease in the uniax- 
ial anisotropy of a freestanding Fe monolayer at a temper- 
ature of about 1000 K. However, the magnetization and 
thus the resulting effective anisotropy will already have 
vanished below 500 K for this case. Therefore the decreas- 
ing magnetization is the main source for the temperature 
dependence of K2(T) in the temperature range we are 
looking at in the present paper. 

In the following, we treat as an example S* = 1, which 
is the lowest spin value with a nontrivial second-order 
anisotropy. For this case, the single-ion decoupling func- 
tion reads ^ — {S'^){S^S^), see equation (||). Then we 
obtain from equation ( |2q ) with n = 1, m = and n = 
l,m = 1, respectively, two coupled equations for (S^) and 
(S^S^), after taking into account the regularity conditions 
of Appendix B for 5' = 1 . The resulting equations are 



(^ 2 jBVZr 
(5^}(l - (Byzf] (2 - iByz)A + 2{S'S') X 



(1 + (Byzy 

^{S'){2-{Byzf^ 



2{3{S^S')-2){l-{B''Z) 



JX r7\2 



^\ + (B-lZf </)(T)=4, (34) 



with 



4>{T) = 



dk^ / dfcj, coth (/3£;k/2) . (35) 



1 + {B^/Zf 



These equations have to be solved numerically in order 
to obtain (5^) and {S^S""). Together with equation @) 
these determine the magnitude M{T) and the orientation 
angle 0o of the magnetization, cf. equations (p9) . 

In the next Section we compare the RPA method with 
a (Bragg- Williams) mean field approximation for the ex- 
change coupling term for a quantum-mechanical spin S = 
1. The resulting expectation values are obtained by diag- 
onalizing the corresponding dynamical matrix consisting 
of the molecular field, the external magnetic field and the 
single- ion anisotropy. For details see Ref. (Q). 



3 Results 

In this section, we display results of our calculations for a 
Heisenberg Hamiltonian plus second-order uniaxial single- 
ion anisotropy for a square monolayer with spin 5' = 1. If 
not stated otherwise, the interactions will be normalized 
to the single-ion anisotropy coefficient at zero temperature 
K2 — K2{T = 0) > 0. We use for the exchange coupling 
J/K2{0) = 100. With zero magnetic field and these pa- 
rameters, the RPA method predicts a Curie temperature 
T^P^/J = 0.989, cf. Appendix A. The mean-field the- 
ory yields with the same parameters a Curie temperature 
j^MFT I J _ 2.667, which is about a factor of three larger. 
The temperature at which the magnetization reaches the 
field direction {{S^) — > 0) is called the reorientation tem- 
perature Tpi, and the corresponding magnetic field the re- 
orientation field Bf^. 

In Fig. 1(a) we show the results obtained from the RPA 
method for the components of the magnetization, {S^} and 
(S^), as functions of the external magnetic field B^ /K2{0) 
in x-direction {B^ = B^ = 0)j for different reduced tem- 
peratures T/Tc- This field- induced magnetic reorienta- 
tion is characterized by a decreasing {S^) and an increas- 
ing {S^). The magnetization reaches the in-plane direc- 
tion {{S^) = 0) at a field strength B^ depending on the 
temperature. For the lowest temperature [T/Tc — 0.21), 
we observe a jump in the components of the magneti- 
zation at the corresponding reorientation field of about 
B^/K2{0) — 0.52. This is probably caused by the par- 
ticular kind of the single-ion anisotropy decoupling ap- 
plied in our procedure. After complete reorientation {S^} 
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Figure 1. Components of the magnetization (S^) (solid lines) 
and {S^) (dashed lines) for different reduced temperatures 
T/Tc as functions of the external magnetic field in x-direction, 
B^ / K2{0) , which is normalized to the anisotropy coefficient at 
temperature T = 0. Results of RPA (a) and MFT (b) calcu- 
lations are compared. The Curie temperatures are T^^ /J = 
0.989 within RPA, and T^^^"^ j J = 2.667 within MFT, using 
J/K2(0) = 100. The reduced temperatures in (a) and (b) are 
chosen in such a way that the magnetizations (S^) are approx- 
imately the same in RPA and MFT at B"" = 0. 



shows a nearly constant behaviour with increasing field. 
In Fig. 1(b) we show the corresponding MFT results with 
the same parameters for J and K2{Q). In order to com- 
pare the different shapes of the magnetization curves in 
RPA and MFT, we have chosen different reduced tem- 
peratures TJTq^'^ for the MFT calculations in such a 
way that the relative magnetizations (S"^) at B^ — are 
about the same in both cases. One observes that (S"^) de- 
creases near B^ more rapidly within the RPA approach 
than within MFT. Also, whereas (5"^) as computed with 
MFT exhibits an almost linear behaviour as a function of 
B^ / K2{0) ^ the corresponding RPA- dependence is clearly 
curved. The main difference is that the reorientation fields 
-B^ are considerably smaller within the RPA approach for 
all temperatures investigated. This we attribute to the 
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Figure 2. The orientation angle 0o of the magnetization is 
shown as a function of B^ / K2{Q) for the same reduced tem- 
peratures as in Fig.l. RPA results (a) are compared with MFT 
results (b). 



fact that a magnetic field has a stronger influence on 2D 
than on 3D systems. Note that the MFT results depend 
only on the coordination number and not on the spatial 
dimensionality, thus MFT handles a 2D system similar to 
a 3D one. 

In Fig.2(a,b) we display the orientation angles 0q{T, 
B^) corresponding to the situation of Fig.l from a per- 
pendicular (6*0 = 0°) to an in- plane direction (6*0 = 90°) 
of the monolayer magnetization. 

The effective anisotropy coefficient, K2{T), is deter- 
mined by equation (B2), which we apply for B^ — 0. 
As mentioned before, this ansatz is physically meaning- 
ful only if the dependence of K2{T) on B^ is small. In 
Fig. 3 the corresponding dependence is shown for different 
temperatures. For all temperatures we obtain only a weak 
dependence of K2{T) on B^ /K2{Q) for small fields, which 
becomes stronger as the reorientation field strength Bj^ 
is approached. Therefore, we have used the small value 
B^/K2{0) = 0.1 to determine K2{T) as a function of 
the reduced temperature. In this case the reorientation 
temperature Tr is close to Tc- The resulting effective 
anisotropy K2{T)/K2{0) and the corresponding orienta- 
tion angle 0o{T) are shown in Fig. 4. K2{T) as obtained 
from RPA is an almost straight line between T = and 
T — Tji ~ Tc- On the other hand, the corresponding 
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Figure 3. Non-perturbative RPA-calculations for the effec- Figure 4. Effective anisotropy coefficient K2{T)/K2{Q) at 

tive single-ion anisotropy coefficient K2{T,B^) normalized to B^/K2{0) = 0.1 and orientation angle Oo{T) are shown as 

K^iO) are shown as a function of the external magnetic field functions of the reduced temperatures T/Tc- Large dots cor- 

B^ /K2{0) for different reduced temperatures T/Tc- respond to RPA, and small squares to MFT results. 



behaviour of K2{T) calculated by MFT is linear at ele- 
vated temperatures but shows an exponential behaviour 
when approaching T — 0. This has the consequence that 
for 5* = 1 the MFT approach yields a considerably smaller 
value of i^2(0) than the RPA method, if observed values of 
K2{T) (measured e.g. at T/Tc ~ 0.7 §|) are extrapolated 
to T = 0. Only at T = are anisotropy constants from ab- 
initio calculations available which can then be compared 
with extrapolated measurements. Note that the range of 
the exponential behaviour of K2{T) when applying MFT 
1^,0 shrinks for large S or for classical spins . 

We do not observe a large difference between RPA and 
MFT results for the orientation angle Oq (T) as a function 
of the reduced temperature T/Tc- Note, however, the dif- 
ferent temperature scale. Tc as obtained from MFT is 2.7 
times larger than the corresponding RPA result for the 
parameters under consideration. 

In Fig. 4 the reorientation temperature is close to the 
Curie temperature, since we have used a small field of 
B^ /K2{0) = 0.1. When applying a larger field, for in- 
stance B^ /K2{Q) — 0.4, the reorientation takes place at 
lower temperatures. This is demonstrated in Fig. 5(a), 
which shows {S^) and {S^) as functions of the reduced 
temperature resulting from RPA calculations, yielding a 
reorientation temperature of about Tji ~ 0.65 T^. The 
corresponding magnitude of the magnetization M{T) = 
y/iS^)"^ + (S*^)^ is also shown, as well as the second mo- 
ments (S^S^) and {S^S'^), which approach the value S'(S'-|- 
l)/3 = 2/3 for large temperatures. For comparison the 
corresponding results of a MFT calculation are displayed 
in Fig. 5(b). For the same applied field the reorientation 
temperature Tji ~ 0.9 Tc is considerably higher in this 
case. We emphasize the long tail in particular of the mag- 
netization M{T) at large temperatures in the RPA calcu- 
lations of Fig. 5(a), which is absent in the MFT results of 
Fig. 5(b). This behaviour is due to the strong effect of ex- 



ternal m agn etic fields on the properties of 2D Heisenberg 
magnets |llO|. 

In Ref. [ p3[ we showed the temperature dependence 
of the single-ion anisotropy coefficients of a ferromagnetic 
monolayer obtained with a thermodynamic perturbation 
theory. The magnetization was calculated in the frame- 
work of RPA, considering the isotropic exchange coupling 
and an external magnetic field B^ in the unperturbed 
Hamiltonian. In perturbation theory, it is unavoidable to 
apply a finite magnetic field in order to obtain a magneti- 
zation at finite temperatures (Mermin- Wagner theorem). 
Thus, the perturbative results for the effective anisotropy 
K2{T) depend sensitively on the strength of the magnetic 
field B"- Even with this approach, we found significant dif- 
ferences from the corresponding mean-field theory calcula- 
tions. As discussed in Section 2, equation ( |32| ) determines 
within RPA the effective single-ion anisotropy K2 (T) non- 
perturbatively. In Fig. 6 we compare the non-perturbative 
result for K2{T), already shown in Fig. 4, with perturba- 
tive results as a function of the reduced temperature. We 
use B^/K2{0) = 0.01, 0.1, and 1.0. As already mentioned 
the absolute value of K2{T) / K2{0) as calculated pertur- 
batively depends considerably on the magnetic field B^ - 
Except for the rounding at elevated temperatures, which 
is clearly an effect of B^ , the shapes of K2 (T) as calculated 
perturbatively and non-perturbatively within the RPA ap- 
proach look similar to each other. The differences between 
perturbative and non-perturbative results become larger 
with increasing -^^2(0). 



4 Discussion and Conclusion 

In the present paper, we have applied many-body Green's 
function theory for the calculation of the magnetic proper- 
ties of a Heisenberg monolayer with second-order uniaxial 
single-ion anisotropy at finite temperatures. This method 
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allows calculations over the entire temperature range of in- 
terest in contrast to other methods, which are only valid 
at low (Holstein-Primakoff approach) or high tempera- 
tures (high temperature expansions). We have used the 
Tyablikov (RPA) decoupling for the exchange interaction 
terms, and the Anderson-Callen decoupling for the anisotr- 
opy terms. For the latter, we have investigated various 
other decoupling schemes, which partly break down at 
higher temperatures or give results similar to the Anderson- 
Callen decoupling, see Appendix A. The results are far 
more sensitive to a variation of the strength of the anisotr- 
opy than to the different decoupling procedures. We em- 
phasize the fact that the present method fulfills the Mer- 
min- Wagner theorem in the limiting case of an isotropic 
2D Heisenberg magnet, in contrast to the mean- field ap- 
proximation applied formerly. 

Our main investigations are concerned with the reori- 
entation of the magnetization induced by a magnetic field 
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perpendicular to the easy axis at finite temperatures, and 
with a non-perturbative calculation of the temperature 
dependent (effective) anisotropy. We have in particular 
investigated the monolayer for spin S ~ 1 which is the 
lowest spin quantum number with a non-trivial second- 
order anisotropy. By solving the equations of motion for 
the Green's functions we calculate the components of the 
magnetization directly, which allows an immediate deter- 
mination of the orientation angle. The effective anisotropy 
coefficient is calculated from the condition that the free 
energy has a minimum at the orientation angle. As dis- 
cussed in Section 2 this is a non-perturbative approach in 
the sense that all quantities entering in the final expression 
are determined by the full Hamiltonian. This is shown to 
be an improvement over the usual thermodynamical per- 
turbative treatment, where the unperturbed part of the 
Hamiltonian must contain the magnetic field. Therefore 
the corresponding results for the effective anisotropy nec- 
essarily depend on the magnetic field. 

Prior to the present work, the magnetic reorientation 
and the effective anisotropics have been calculated mainly 
within the framework of mean-field theory (MFT) |§|,0|. 
The magnetic reorientation has also been investigated with 
other theoretical methods ||l^,^. In the presence of aniso- 
tropic interactions, MFT is expected to yield qualitatively 
correct results for 2D magnetic systems. By comparison, 
however, we find significant quantitative differences be- 
tween the results obtained with MFT and the present 
Green's function theory. Owing to the magnetic fluctua- 
tions, one finds a different temperature scale with RPA as 
compared to MFT. With the same input parameters, the 
Curie temperature Tc of the present 2D system is about a 
factor of three larger in MFT than in RPA. If the temper- 
ature is rescaled with respect to the corresponding Curie 
temperatures the orientation angles of the magnetization 
look very similar in RPA and MFT. This we consider to be 
a non-trivial result because one would not expect this from 
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the very beginning. The effective anisotropies, however, 
behave differently, also when rescaling the temperature. 
When using RPA the temperature dependence of the effec- 
tive anisotropy coefficient behaves linearly over the whole 
temperature range < T < Tc- fn MFT one observes an 
exponential dependence at low temperatures and small 
S. This different behaviour would manifest itself when 
extrapolating measurements performed at finite tempera- 
tures down to r = in order to compare with ab-initio 
calculations [B for the anisotropy coefficient K2{T — 0), 
which are available only there. 

One also finds in RPA a stronger influence of the ex- 
ternal magnetic field on the reorientation of the magne- 
tization, which is reminiscent of 2D Heisenberg magnets 
]fO| . We observe that, at a fixed temperature, a weaker 
reorientation field B^ is required in RPA than in MFT in 
order to align the magnetization along the field direction. 
At a fixed magnetic field, on the other hand, one obtains 
in RPA a lower reduced reorientation temperature T^jTc- 
Furthermore, the magnetization calculated by RPA has a 
long tail at large temperatures due to the magnetic field, 
which is absent in MFT results. In general, one expects 
the differences between RPA and MFT to decrease as the 
number of film layers increases. 

Investigations into extending the model are in progress. 
The other decoupling procedure different from that of 
Anderson-Callen, which also works up to Tc-, see Fig. 7, 
will be investigated. In this case, none of the eigenval- 
ues of the secular problem vanishes and one has to mod- 
ify the formalism appropriately. A next step is to include 
the magnetic dipole coupling competing with the uniax- 
ial single-ion anisotropy. The effect of the dipole coupling 
will become more important when the calculations are ex- 
tended to magnetic films with several layers in order to 
treat also the dependence of the reorientation on the film 
thickness. In this case, the layer dependent anisotropies 
might lead to a temperature driven reorientation with- 
out the application of an external magnetic field, since 
the temperature dependence of the anisotropies will be 
different for surface and interior film layers. In addition, 
calculations for larger spins will be attempted, at least for 
spin 5 = 2, in order to be able to treat the fourth-order 
single-ion anisotropy K/^. Then one can set up a phase dia- 
gram e.g. in the Ki — iir4-plane, which will show the region 
of stable magnetization directions and the location of the 
temperature driven magnetic reorientation. Furthermore, 
one can investigate whether the magnetic reorientation 
takes place continuously or discontinuously |^,0. 



Appendix A: Decoupling schemes for the single- 
ion anisotropy term 

In this appendix we discuss various decoupling procedures 
for the single-ion anisotropy term, and give arguments for 
using the Anderson-Callen decoupling. For an anisotropy 
strong compared to the exchange coupling, a method us- 
ing different Green's functions has been applied formerly 
for spin 5* = I |24]. However, this method leads to an 



overdetermined system of equations for the expectation 
values (5^) and (S^S^). On the other hand, for small 
anisotropies, as is the case for the magnetic systems con- 
sidered in the present paper, different approaches have 
been proposed by Anderson-Callen [^ and by Lines [ p2[ , 
which we shall treat here. 

For the decoupling of the higher-order Green's func- 
tions coming from the single-ion anisotropy term, an ansatz 
of the following form is used: 



{{StSf + StS^;{S]r{S-r)) 



^±,mn ^±mn ^ ^gg^ 



where the functions ^^ '™" have to be determined. For 
m — 0, n = 1, Lines [g2| proved that <P^ can be ex- 
pressed by the following ratios of commutators 



^: 



,01 



{[Sr,S+{2Sf + 1)]) _ 3((5f )2) - S{S + 1) 



For m = n = 1 one obtains 



(S!) 



<^; 



_,ii ^ iSS!Sr,S+(2S! + l)\) 



(37) 



(38) 



8((gf )3) - 3{{Sfr) + (1 ~ 4^(5 + l)){Sf) + S{S + 1) 

Note that the t^v™" depend only on a single lattice site 
i, as expected. Similar expressions may be derived for 

An alternative expression for <P^ '™"' can be obtained 
by replacing the Green's functions in equation ( pq ) by 
their respective expectation values, resulting in 

_ (39) 

One obtains for n = 1 equation ( p^ ) shown below. For 
spin 5 = 1, which we treat in the present paper, we need 
m = 0, 1 and n = 1. One has S{S + 1) = 2, {S^f = Sf , 
and (5*1 )■* = (^f )^, thus the ratios of commutators are 



<?: 



+,01 ^ 3((5f)^)-2 



<P 



-i-,ii 



-1, 



(41) 



and the ratios of the expectation values are given by 
2 + {St) - 3((5f )2) 



<p: 



{st)-{{sm ' 



<pi 



-1. 



(42) 



Note that Lines [g2| mixes both approaches, since he uses 
equation ( |3^ ) for even m and equation (40) for odd m. 

In order to investigate the differences between these 
decoupling procedures, we treat the special case B^ = 
By = 0, for which {Sf) = (Sf) = 0. Then one has to 
determine the expectation values (Sf) and (SfSf). For 
S — 1 we obtain from equations ( |33| ) and (p4), dropping 
the site index i for the ferromagnetic system, 






(5^)(I + 20o) 
1 + 200 



1 + 2(j)i + (j)o + 3(j)i(j)o 



(43) 
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S{S + l)((Sf)™> - ((5'f)'"+i> - ((S-f )-+2> 
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Figure 7. The magnetization (S^) and the second moment 
{S^S'^) are shown as function of the temperature for various 
decouplin g s chemes discussed in the text, (a) decoupling of 
equation ( |4l| ) ; ( b) decoupling of equation (tt2|); (c) decoupling 
of equation (M); (d) Anderson-Callen decoupling ( equation 
(Eq), small dots connected by thin line). Only the decouplings 
(cj and (d) work up to {S^) = 0, i.e. up to the Curie tempera- 
ture Tc, whereas (a) and (b) fail for T < Tc- We include also 
a result with Anderson-Callen decoupling and a magnetic field 
ByK2{0) = 1. 



with 



= ^J2(^MPE^)-l 



E^ = B' + {S')J{A - 7k) + K2 ^" 



(44) 
(45) 



with m = 0, 1, where <?™ is either c^+'^i or S+'"^^. 

When solving; these equations for the decouphng pro- 
cedures (fill), ([42|), we find that both procedures fail at 
temperatures far below Tc, cf. Fig. 7 (situation (a),(b)). 
Thus one cannot use both decouplings up to the Curie 
temperature Tc- 

When calculating Curie temperatures, Lines ^^ has 
circumvented the problem associated with equations ( [43[ ) 
by (inconsistently) using in the dispersion relations E^^, 
equation (^5|), expectation values {S^) and {S'^S^) as ob- 
tained from the one-particle density operator, using a the- 
orem originating from Callen and Strikman [ P5| . 

If, however, one uses 



<?+ 



,ml 



(?)+.oi Cor ^^ 



, r??. 1 



c|>+ ■' \ respectively ) (46 ) 



for m = 0, 1, the procedure works well and yields rea- 
sonably results for {S^) and {S^S^) between T = and 
T — Tc, cf. Fig. 7 (situation (c)), and the following discus- 
sion about the Curie temperature Tc- 



In this case, (/iq = 0i = , see equation (Q), and 
one finds that the procedures for t^+'^i and <?i+'"^ give the 
same result. This can be understood by equating ^+'°i = 
^+^01^ which yields the second of equations ( p3| ) determin- 
ing(5^^). 

Using this decoupling in the Green's function formal- 
ism one can derive expressions for the ^~'™", using also 
^-,01 _ (^-41 xhis leads to a 3 x 3 secular problem for 
equation dq) with three non-vanishing eigenvalues, since 
here ^+'™" ^ <p-,mn_ ^^ ^-j-^^g ^^^gg |-j^g regularity condition 
( pi|) cannot be applied. This case is more complicated to 
handle than the Anderson-Callen decoupling, in which one 
of the eigenvalues vanishes, cf. equation (p4). 

The Anderson-Callen decoupling is based on the paper 
by Callen |23| , in which correlations beyond the RPA are 
included in the decoupling of the Green's functions of the 
exchange terms. In this case, the essential Green's function 
is non-diagonal (i ^ j) and the decoupling reads (e.g. for 
m — 0, n — I) 

{{StSf;S-)) ^ {S^){{sf;S-))-^-^{Srsf){{St;S^)) . 

(47) 
Neglecting the second term corresponds to the RPA de- 
coupling. 

The proposal of Anderson and Callen [O is to use the 
same decoupling also for Green's functions with i — j, 
resulting from the anisotropy terms 






'^^{{sr-^s-}) 



(48) 



In this case, one does not distinguish between different 
m. Since here ^+ ~ (P^ ~ (P, which is the quantity <!> 
given in equation (|9|) , one of the eigenvalues turns out to 
be zero, cf. equation (M). This is a prerequisite for being 
able to use the procedure outlined in Section 2 for the 
determination of the expectation values {S^) and {S^S^}. 

Inspecting Fig. 7 shows that the two decouphngs, 
which work up to the Curie temperature, do not give very 
different results. Therefore, we adopt the decoupling which 
is easier to apply: the Anderson-Callen decoupling. 

In the remainder of this appendix, we derive expres- 
sions which determine the Curie temperature, Tc{J, K2), 
for the case 5 = 1 and (l>+'°^ = c^+'H = c^, cf. equa- 
tion (^). The consideration of different c^+'O^ ^ <P+'^^, 
cf. equation ([4l|), leads to severe complications whilst de- 
termining Tc, or, as we have seen, when calculating the 
magnetization (S*^) as function of T. 

Tc is calculated from performing the limit (S*^) -^ 
in equations (|4^). To lowest order in (S^), one obtains for 
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Figure 8. Curie temperatures normalized to the exchange 
coupling, Tc / J, are shown as a function of the anisotropy cou- 
pling strength K2{0)/J, using equation (BSl) with a — 3/4 
(decoupling according to equation (jAw) and a = 2/3 corre- 
sponding to the Anderson- C alien decoupling (equation (f4q) . 
The Curie temperatures at 7^2(0)/ J = 0.01 correspond to 
those of Fig.7 (cases (c) and (d)). Also, the results from the ap- 
proximate expression, equation (M), for a = 3/4 and a — 2/3 
are included in the figure. 



T^Tc 

{{S'f){T < Tc)- 2/3 +(S')y 4, 
and for the decoupling functions 

•PiTc) = a (S^) , 



(49) 



(50) 



with a = 3/4 for the decoupling given by equation (Bq) 
and a — 2/3 for the Anderson-Callen decoupling, equa- 
tions (^ or (^ ) . Close to Tc one also obtains from equa- 
tion (kH) 



0(Tc) 



Tc 

{S') N 



^^(j(4-27k) + i^2> 



and from equation ( p3| ) 

HTc) 



3(5^ 



(51) 



(52) 



Combining these equations and converting the sum into 
an integral gives 



Tc = 



8^2 






dk„ 



1 



^(4 - 7k) + K2a 



(53) 

The Curie temperatures calculated with this formula are 
displayed in Fig. 8 for the two single-ion decouplings under 
consideration. 

Since this integral is dominated by small wave numbers 
it can be approximately evaluated when expanding 7k up 
to the leading order in k: 



Tc 



?7rJ/3 



In(l + 27r2j/if2a) 



(54) 
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Figure 9. The temperature dependence of (S^) (a) and 

(S^S'^) (b) for different strengths of the anisotropy coefficient 
K2{0)/J = 0.001, 0.01, 0.05, 0.1, 0.2, and 0.5 is displayed. 



This expression can be used as a quick estimate for the 
Curie temperature. It overestimates the result as obtained 
from equation (psh by less than 10%. The corresponding 
results are also shown in Fig. 8. 

We emphasize the fact that the various decoupling pro- 
cedures do not yield very different results for the magne- 
tization, see Fig.7, or the Curie temperature Tc, cf. Fig. 8. 
We note that the results are much more sensitive to vary- 
ing the strength of the anisotropy coefficient K2 (0) than 
to the different decoupling procedures. This is shown in 
Fig. 9 by plotting the magnetization (S"^) and (S^S^) as 
function of the temperature for different values K2{0)/J 
of the anisotropy coefficient normalized to the exchange 
coupling strength J for the case of the Anderson-Callen 
decoupling. 



Appendix B 

appb) In this Appendix we list explicitly a number of rela- 
tions obtained from the regularity condition equation (12 



2ZA': 



^A 



-\-,7nn Tj — 



a: 



1^4 



(54) 



Remember that we have used the notations B^ — B^ ± 
iBy,Z = B'+K2'P, and AX" = ([S*", (S'^)™(S-)"]_i). 
In Table ^ we tabulate some commutator relations. 
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m 


n 


Al'i'"" 


^+,rnn 


Azr" 





1 


~{S~) 


2{Sn 










2 


-2{S-S~) 


4{S"S-) + 2(S-) 







1 


1 


~{s^s-) 


3{5"S") - (S"> - 5(5 + 1) 


(S-5-) 







3 


-3{S-S-S-} 


6{S''S~S-)+6{S-S-) 







1 


2 


~2{S^S-S-) 


5{S"S"S-> + (S^S-) - {S^)S{S + 1) 


(5-5-5- 


) 


2 


1 


-(S"S"S-) 


+(1 - 25(5- + 1)){S"> + 5(5 + 1) 


2(5^5-5 


-) + {s-s-) 



Table 1. Inhomogeneities A^'^" = (f5", (5^)™(5-)"]-i), a = +,-,2, of the anticommu- 
tator Green's functions, of. equation (Q), for different integers m and n. 



(5" 5- 5-) = 



B-Z 



2^2 - 3B+B 

-('3{5"5")- 5(5+1)^ 
(B-f 



4(5"5"5"> - (5"> ('25(5 + 1) - 1) 
2^2-3B+B- B+B 



2^2 _ B+B- 



^2 



{S''S^S^} + ^{S^) 



2^2 - 3B+B- 
Z(2^2 -SB+B- 



5{5"5"5"> - (5"> 35(5 + 1) - 1 - 3(5" 5"> - 3(3 + 1) 



2Z'^ - 3B+B' 



2^2 - B+B- 



j(5"5"5")- (5") (35(5 + 1) 



(57) 



We obtain with m = 0,n = I 

with m — 0,n — 2 and 777, = 1, n = 1 
(i3-)2 



(5-5-) = 
(5^5-) = 



(55) 



-(3(5^5^) -5(5 + 1) 
(3(5^5") -5(5 + 1)) 
(5^) , (56) 



2Z2 - B+B 

B-Z 
2Z^~~B+B 

B- 
~2Z 



and with m — 0,n — 3, and ?77 = 1, n = 2, and m = 2,n 
1 given in equations (B^ displayed above. 
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